A world record in Atlantic City 
and the length of the shooter's hand at craps 
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It was widely reported in the media that, on 23 
May 2009, at the Borgata Hotel Casino & Spa in At- 
lantic City, Patricia DeMaurcQ, playing craps for only 
the second time, rolled the dice for four hours and 18 
minutes, finally scvening out at the 154th roll. Ini- 
tial estimates of the probability of this event ranged 
from one chance in 3.5 billion [5] to one chance in 1.56 
trillion [llj . Subsequent computations agreed on one 
chance in 5.6 (or 5.59) billion [2H3QJJ. 

This established a new world record, previously 
held by the late Stanley Fujitakc (118 rolls, 28 May 
1989, California Hotel and Casino, Las Vegas) [T]. 
One might ask how reliable these numbers (118 and 
154) are. In Mr. Fujitake's case, casino personnel re- 
played the surveillance videotape to confirm the num- 
ber of rolls and the duration of time (three hours and 
six minutes) . We imagine that the same happened in 
Ms. DeMauro's case. 

There is also a report that Mr. Fujitake's record 
was broken earlier by a gentleman known only as 
The Captain (148 rolls, July 2005, Atlantic City) 
[HI Part 4]. However, this incident is not well doc- 
umented and was unknown to Borgata officials. In 
fact, a statistical argument has been offered [H p. 480] 
suggesting that the story is apocryphal. 

Our aim in this article is not simply to derive a 
more accurate probability, but to show that this ap- 
parently prosaic problem involves some interesting 
mathematics, including Markov chains, matrix the- 
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ory, and Galois theory. 

Background 

Craps is played by rolling a pair of dice repeatedly. 
For most bets, only the sum of the numbers appearing 
on the two dice matters, and this sum has distribution 

j ^ 1 j = 2, 3,..., 12. (1) 

The basic bet at craps is the pass-line bet, which is 
defined as follows. The first roll is the come-out roll. 
If 7 or 11 appears (a natural), the bettor wins. If 2, 
3, or 12 appears (a craps number), the bettor loses. 
If a number belonging to 

& := {4,5,6,8,9,10} 

appears, that number becomes the point. The dice 
continue to be rolled until the point is repeated (or 
made), in which case the bettor wins, or 7 appears, in 
which case the bettor loses. The latter event is called 
a seven out. A win pays even money. The first roll 
following a decision is a new come-out roll, beginning 
the process again. 

A shooter is permitted to roll the dice until he or 
she sevens out. The sequence of rolls by the shooter 
is called the shooter's hand. Notice that the shooter's 
hand can contain winning 7s and losing decisions 
prior to the seven out. The length of the shooter's 
hand (i.e., the number of rolls) is a random variable 
we will denote by L. Our concern here is with 

t(ri) := P{L > n), n > 1, (2) 

the tail of the distribution of L. For example, £(154) 
is the probability of achieving a hand at least as long 
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as that of Ms. DeMauro. As can be easily verified 
from ©, ©, or © below, t(154) « 0.178 882 426 x 
10 -9 ; to state it in the way preferred by the media, 
this amounts to one chance in 5.59 billion, approxi- 
mately The 1 in 3.5 billion figure came from a simu- 
lation that was not long enough. The 1 in 1.56 trillion 
figure came from (1 — 7T7) 154 , which is the right answer 
to the wrong question. 



Two methods 

We know of two methods for evaluating the tail prob- 
abilities ([2|) . The first is by recursion. As pointed out 
in [3], t(l) = t(2) = 1 and 



matrix, which can be inferred from (TT]), is 



t(n) 



1 - ^jt(n - 1) + 5>i(l - - n 7 ) n -" 

n-1 

+ E *i X> - *i - ^y- 2 ^ t(n - I) (3) 
je&> 1=2 



for each n > 3. Indeed, for the event that the shooter 
sevens out in no fewer than n rolls to occur, consider 
the result of the initial come-out roll. If a natural or 
a craps number occurs, then, beginning with the next 
roll, the shooter must seven out in no fewer than n— 1 
rolls. If a point number occurs, then there are two 
possibilities. Either the point is still unresolved after 
77,-2 additional rolls, or it is made at roll I for some 
I € {2,3, ...,77— 1} and the shooter subsequently 
sevens out in no fewer than 77 — I rolls. 

The second method, first suggested, to the best of 
our knowledge, by Peter A. Griffin in 1987 (unpub- 
lished) and rediscovered several times since, is based 
on a Markov chain. The state space is 

S := {co, p4-10, p5-9, p6-8, 7o} = {1, 2, 3, 4, 5}, (4) 

whose five states represent the events that the shooter 
is coming out, has established the point 4 or 10, has 
established the point 5 or 9, has established the point 
6 or 8, and has sevened out. The one-step transition 
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/12 6 8 10 \ 

3 27 6 

4 26 6 

5 25 6 
\ 36/ 



(5) 



The probability of sevening out in 77 — 1 rolls or less 
is then just the probability that absorption in state 
7o occurs by the (n— l)th step of the Markov chain, 
starting in state co. A marginal simplification results 
by considering the 4 by 4 principal submatrix Q of 
([5]) corresponding to the transient states. Thus, we 
have 



t(n) 



)1,5 



- 1 ) 



(6) 



Clearly, ([3]) is not a closed-form expression, and we 
do not regard ([6]) as being in closed form either. Is 
there a closed-form expression for t(n)7 

Positivity of the eigenvalues 

We begin by showing that the eigenvalues of Q arc 
positive. The determinant of 



Q-zJ 
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/l2-36z 6 8 10 

3 27- 36z 

4 26 - 36z 
\ 5 25 - 36z/ 



is unaltered by row operations. From the first row 
subtract 6/ (27 — 36z) times the second row, 8/(26 — 
36z) times the third row, and 10/(25 — 36z) times 
the fourth row, cancelling the entries 6/36, 8/36, and 
10/36 and making the (1,1) entry equal to 1/36 times 



12-36z-3- 
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26 - 36z 



25 - 36z 
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The determinant of Q — zl, and therefore the char- 
acteristic polynomial q(z) of Q is then just the prod- 
uct of the diagonal entries in the transformed matrix, 
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which is multiplied by (27 - 36z)(26 - 36z)(25 - 
36z)/(36) 4 . Thus, 

q(z) = [(12 - 36z)(27 - 36z)(26 - 36z)(25 - 36z) 

- 18(26 - 36z)(25- 36z) 

- 32(27- 36z)(25- 36z) 

- 50(27 - 36z)(26 - 36z)]/(36) 4 . 

Wc find that g(l),g(27/36),g(26/36),g(25/36),g(0) 
alternate signs and therefore the eigenvalues are pos- 
itive and interlaced between the diagonal entries (ig- 
noring the entry 12/36). More precisely, denoting the 
eigenvalues by 1 > e\ > e 2 > > e 4 > 0, we have 

27 26 25 

1>ei> 36 >e2> 36 >e3> 36 >e4>0 ' 

The matrix Q, which has the structure of an ar- 
rowhead matrix [5], is positive definite, although 
not symmetric. This is easily seen by applying the 
same type of row operations to the symmetric part 
■A = 7j(Q + Q J ) of Q to show that the eigenvalues of 
A interlace its diagonal elements (except 12/36). But 
a symmetric matrix is positive definite if and only if 
all its eigenvalues are positive, and a non-symmetric 
matrix is positive definite if and only if its symmetric 
part is positive definite, confirming that Q is positive 
definite. 



A closed-form expression 

The eigenvalues of Q arc the four roots of the quartic 
equation q(z) = or 

23328z 4 - 58320z 3 + 51534z 2 - 18321z + 1975 = 0, 

while P has an additional eigenvalue, 1, the spectral 
radius. We can use the quartic formula (or Mathe- 
matica) to find these roots. We notice that the com- 
plex number 

a]=(1/3 9829 



(1/3 



appears three times in each root. Fortunately, a is 
positive, as we see by writing ( in polar form, that 
is, £ = re 10 . Wc obtain 



2V9829 cos 
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The four eigenvalues of Q can be expressed as 

e(l,l), 
e(l,-l), 
e(-l,l), 
e(-l,-l), 



ei 

(•■2 
f'4 



where 

e(u, v) 




2136u 



349 + a 



Next we need to find right eigenvectors correspond- 
ing to the five eigenvalues of P. Fortunately, these 
eigenvectors can be expressed in terms of the eigen- 
values. Indeed, with r(x) defined to be the vector- 
valued function 

/ -5 + (l/5)a: \ 

-175 + (581/15)x - (21/10)a; 2 + (l/30)a; 3 
275/2 - (1199/40):r + (8/5).t 2 - (l/40)x 3 
1 

V o / 

we find that right eigenvectors corresponding to 
eigenvalues 1, ei, e%, e^, e 4 are 

(1,1,1,1,1) T , r(36 ei ), r(36e 2 ), r(36e 3 ), r(36e 4 ), 

respectively. Letting R denote the matrix whose 
columns are these right eigenvectors and putting 
L := it -1 , the rows of which are left eigenvectors, 
we know by ([6|) and the spectral representation that 



t{n)^l-{Rdmg{l,e'l- 1 ,el 



-1 n-l 

' 3 ' e 



1,5- 



where 



After much algebra (and with some help from 
Mathematica), we obtain 



C := -710369 + 18iVl373296647, 



t(n) = c x e r l 1 +c 2 e£ 1 + c 3 e£" 



c 4 e 4 



(8) 
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where the coefficients are defined in terms of the 
eigenvalues and the function 

f(w, x, y, z) := (-25 + 36w)[4835 - 5580(x + y + z) 
+ 6480{xy + xz + yz) - 7776xyz] 
/[38880(«; - a:) (to - y)(w - z)\ 



as follows: 



Cl 
f'2 
C3 
'■'4 



/(ei,e 2 ,e3,e 4 ), 
/(e2,e 3 ,e 4 ,ei), 

/(es, 64, 61,62), 

/(e 4 ,ei,e 2 , e 3 ). 



Of course, ijSJ) is our closed-form expression. 

Incidentally, the fact that t(l) = i(2) = 1 implies 
that 

ci + c 2 + C3 + c 4 = 1 (9) 

and 

ciei + c 2 e 2 + c 3 e 3 + c 4 e 4 = 1. 

In a sequence of independent Bernoulli trials, each 
with success probability p, the number of trials X 
needed to achieve the first success has the geometric 
distribution with parameter p, and 



P(X > n) = (l-p) 



n-l 



n > 1. 



It follows that the distribution of L is a linear com- 
bination of four geometric distributions. It is not a 
convex combination: j9]) holds but, as we will see, 

Ci > 0, c 2 < 0, c 3 < 0, c 4 < 0. 

In particular, we have the inequality 



t(n) < cie" 



n-l 



n > 1, 



as well as the asymptotic formula 

t(n) ~ cie" -1 as n — > 00. 



(10) 
(11) 



Numerical approximations 

Rounding to 18 decimal places, the nonunit eigenval- 
ues of P arc 

ei « 0.862 473 751 659 322 030, 



e 2 ps 0.741 708 271 459 795 977, 
e 3 ps 0.709 206 775 794 379 015, 
e 4 ps 0.186 611 201 086 502 979, 



and the coefficients in 



arc 



ci ps 1.211 844 812 464 518 572, 
c 2 w -0.006 375 542 263 784 777, 
c 3 ss -0.004 042 671 248 651 503, 
c 4 ss -0.201426 598 952 082 292. 

These numbers will give very accurate results over a 
wide range of values of n. 

The result (fTTj) shows that the leading term in © 
may be adequate for large n; it can be shown that 

1 < de? _1 /*(n) < 1 + 10~ m 

for m = 3 if n > 19; for m — 6 if n > 59; for m = 9 
if n > 104; and for m = 12 if n > 150. 



Crapless craps 

In crapless craps [5J p. 354], as the name suggests, 
there are no craps numbers and 7 is the only natural. 
Therefore, the set of possible point numbers is 

^0 := {2,3,4,5,6,8,9,10,11,12} 

but otherwise the rules of craps apply. 

With Lq denoting the length of the shooter's hand, 
the analogues of (J!])-© are 

So := {co,p2-12,p3-ll,p4-10,p5-9,p6-8,7o} 
= {1,2,3,4,5,6,7}, 



/6 2 4 6 8 10 \ 

1 29 6 

2 28 6 

3 27 6 

4 26 6 

5 25 6 
\0 36/ 
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and 



i (n) :=P(L > n) = 1 - (P^hj- 
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There is an interesting distinction between this 
game and regular craps. The nonunit eigenvalues of 
Po are the roots of the sextic equation 

15116544z 6 - 59206464z 5 + 93137040z 4 
- 73915740z 3 + 30008394z 2 - 5305446z + 172975 
= 0, 

and the corresponding Galois group is, according to 
Maple, the symmetric group S"g. This means that our 
sextic is not solvable by radicals. Thus, it appears 
that there is no closed-form expression for to(n). 

Nevertheless, the analogue of ([5]) holds (with six 
terms). All nonunit eigenvalues belong to (0,1) and 
all coefficients except the leading one are negative. 
Thus, the analogues of (fTQ|) and (fTTj) hold as well. 
Also, the distribution of Lq is a linear combination 
of six geometric distributions. These results are left 
as exercises for the interested reader. 

Finally, i (154) w 0.296 360 068 x 10" 10 , which is 
to say that a hand of length 154 or more is only about 
one-sixth as likely as at regular craps (one chance in 
33.7 billion, approximately). 
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